Statistical Inference of Selection and Divergence from a Time-Dependent Poisson Random Field Model

We apply a recently developed time-dependent Poisson random field model to aligned DNA sequences from two related biological species to estimate selection coefficients and divergence time. We use Markov chain Monte Carlo methods to estimate species divergence time and selection coefficients for each locus. The model assumes that the selective effects of non-synonymous mutations are normally distributed across genetic loci but constant within loci, and synonymous mutations are selectively neutral. In contrast with previous models, we do not assume that the individual species are at population equilibrium after divergence. Using a data set of 91 genes in two Drosophila species, D. melanogaster and D. simulans, we estimate the species divergence time (or 1.68 million years, assuming the haploid effective population size years) and a mean selection coefficient per generation . Although the average selection coefficient is positive, the magnitude of the selection is quite small. Results from numerical simulations are also presented as an accuracy check for the time-dependent model

No severe and global X chromosome inactivation in meiotic male germline of Drosophila

This article is a response to Vibranovski et al

Modularity and Intrinsic Evolvability of Hsp90-Buffered Change

Hsp90 controls dramatic phenotypic transitions in a wide array of morphological features of many organisms. The genetic-background dependence of specific abnormalities and their response to laboratory selection suggested Hsp90 could be an ‘evolutionary capacitor’, allowing developmental variation to accumulate as neutral alleles under normal conditions and manifest selectable morphological differences during environmental stress. The relevance of Hsp90-buffered variation for evolution has been most often challenged by the idea that large morphological changes controlled by Hsp90 are unconditionally deleterious. To address this issue, we tested an Hsp90-buffered abnormality in Drosophila for unselected pleiotropic effects and correlated fitness costs. Up to 120-fold differences in penetrance among six highly related selection lines, started from an initially small number of flies and rapidly selected for and against a deformed eye trait (dfe), did not translate into measurable differences in any of several tests of viability, lifespan or competitive fitness. Nor were 17 dfe Quantitative Trait Loci (QTL) associated with fitness effects in over 1,400 recombinant lines. Our ability to detect measurable effects of inbreeding, media environment and the white mutation in the selection line backgrounds independent of dfe penetrance suggests that, within the limitations of laboratory tests of fitness, this large morphological change controlled by Hsp90 was selectable independent of strong, correlated and unconditionally deleterious effects—abundant, polygenic variation hidden by Hsp90 allows potentially deleterious alleles to be readily replaced during selection by less deleterious alleles with similar phenotypic effects. Hsp90 links environmental stress with the expression of developmental variation controlling unprecedented morphological plasticity. As outlined here and in the companion paper of this issue, the complex genetic architecture of Hsp90-buffered variation supports a remarkable modularity of Hsp90 effects on quantitative and qualitative phenotypes, consistent with the ‘Hsp90 capacitor hypothesis’ and standard quantitative genetic models of threshold traits

Sex-specific Trans-regulatory Variation on the Drosophila melanogaster X Chromosome

The X chromosome constitutes a unique genomic environment because it is present in one copy in males, but two copies in females. This simple fact has motivated several theoretical predictions with respect to how standing genetic variation on the X chromosome should differ from the autosomes. Unmasked expression of deleterious mutations in males and a lower census size are expected to reduce variation, while allelic variants with sexually antagonistic effects, and potentially those with a sex-specific effect, could accumulate on the X chromosome and contribute to increased genetic variation. In addition, incomplete dosage compensation of the X chromosome could potentially dampen the male-specific effects of random mutations, and promote the accumulation of X-linked alleles with sexually dimorphic phenotypic effects. Here we test both the amount and the type of genetic variation on the X chromosome within a population of Drosophila melanogaster, by comparing the proportion of X linked and autosomal trans-regulatory SNPs with a sexually concordant and discordant effect on gene expression. We find that the X chromosome is depleted for SNPs with a sexually concordant effect, but hosts comparatively more SNPs with a sexually discordant effect. Interestingly, the contrasting results for SNPs with sexually concordant and discordant effects are driven by SNPs with a larger influence on expression in females than expression in males. Furthermore, the distribution of these SNPs is shifted towards regions where dosage compensation is predicted to be less complete. These results suggest that intrinsic properties of dosage compensation influence either the accumulation of different types of trans-factors and/or their propensity to accumulate mutations. Our findings document a potential mechanistic basis for sex-specific genetic variation, and identify the X as a reservoir for sexually dimorphic phenotypic variation. These results have general implications for X chromosome evolution, as well as the genetic basis of sex-specific evolutionary change

Mechanisms of gap gene expression canalization in the Drosophila blastoderm

<p>Abstract</p> <p>Background</p> <p>Extensive variation in early gap gene expression in the <it>Drosophila </it>blastoderm is reduced over time because of gap gene cross regulation. This phenomenon is a manifestation of canalization, the ability of an organism to produce a consistent phenotype despite variations in genotype or environment. The canalization of gap gene expression can be understood as arising from the actions of attractors in the gap gene dynamical system.</p> <p>Results</p> <p>In order to better understand the processes of developmental robustness and canalization in the early <it>Drosophila </it>embryo, we investigated the dynamical effects of varying spatial profiles of Bicoid protein concentration on the formation of the expression border of the gap gene <it>hunchback</it>. At several positions on the anterior-posterior axis of the embryo, we analyzed attractors and their basins of attraction in a dynamical model describing expression of four gap genes with the Bicoid concentration profile accounted as a given input in the model equations. This model was tested against a family of Bicoid gradients obtained from individual embryos. These gradients were normalized by two independent methods, which are based on distinct biological hypotheses and provide different magnitudes for Bicoid spatial variability. We showed how the border formation is dictated by the biological initial conditions (the concentration gradient of maternal Hunchback protein) being attracted to specific attracting sets in a local vicinity of the border. Different types of these attracting sets (point attractors or one dimensional attracting manifolds) define several possible mechanisms of border formation. The <it>hunchback </it>border formation is associated with intersection of the spatial gradient of the maternal Hunchback protein and a boundary between the attraction basins of two different point attractors. We demonstrated how the positional variability for <it>hunchback </it>is related to the corresponding variability of the basin boundaries. The observed reduction in variability of the <it>hunchback </it>gene expression can be accounted for by specific geometrical properties of the basin boundaries.</p> <p>Conclusion</p> <p>We clarified the mechanisms of gap gene expression canalization in early <it>Drosophila </it>embryos. These mechanisms were specified in the case of <it>hunchback </it>in well defined terms of the dynamical system theory.</p

Networked buffering: a basic mechanism for distributed robustness in complex adaptive systems

A generic mechanism - networked buffering - is proposed for the generation of robust traits in complex systems. It requires two basic conditions to be satisfied: 1) agents are versatile enough to perform more than one single functional role within a system and 2) agents are degenerate, i.e. there exists partial overlap in the functional capabilities of agents. Given these prerequisites, degenerate systems can readily produce a distributed systemic response to local perturbations. Reciprocally, excess resources related to a single function can indirectly support multiple unrelated functions within a degenerate system. In models of genome:proteome mappings for which localized decision-making and modularity of genetic functions are assumed, we verify that such distributed compensatory effects cause enhanced robustness of system traits. The conditions needed for networked buffering to occur are neither demanding nor rare, supporting the conjecture that degeneracy may fundamentally underpin distributed robustness within several biotic and abiotic systems. For instance, networked buffering offers new insights into systems engineering and planning activities that occur under high uncertainty. It may also help explain recent developments in understanding the origins of resilience within complex ecosystems. \ud \u

Chronic, low-dose rotenone reproduces Lewy neurites found in early stages of Parkinson's disease, reduces mitochondrial movement and slowly kills differentiated SH-SY5Y neural cells

<p>Abstract</p> <p>Background</p> <p>Parkinson's disease, the most common adult neurodegenerative movement disorder, demonstrates a brain-wide pathology that begins pre-clinically with alpha-synuclein aggregates ("Lewy neurites") in processes of gut enteric and vagal motor neurons. Rostral progression into substantia nigra with death of dopamine neurons produces the motor impairment phenotype that yields a clinical diagnosis. The vast majority of Parkinson's disease occurs sporadically, and current models of sporadic Parkinson's disease (sPD) can utilize directly infused or systemic neurotoxins.</p> <p>Results</p> <p>We developed a differentiation protocol for human SH-SY5Y neuroblastoma that yielded non-dividing dopaminergic neural cells with long processes that we then exposed to 50 nM rotenone, a complex I inhibitor used in Parkinson's disease models. After 21 days of rotenone, ~60% of cells died. Their processes retracted and accumulated ASYN-(+) and UB-(+) aggregates that blocked organelle transport. Mitochondrial movement velocities were reduced by 8 days of rotenone and continued to decline over time. No cytoplasmic inclusions resembling Lewy bodies were observed. Gene microarray analyses showed that the majority of genes were under-expressed. qPCR analyses of 11 mtDNA-encoded and 10 nDNA-encoded mitochondrial electron transport chain RNAs' relative expressions revealed small increases in mtDNA-encoded genes and lesser regulation of nDNA-encoded ETC genes.</p> <p>Conclusion</p> <p>Subacute rotenone treatment of differentiated SH-SY5Y neuroblastoma cells causes process retraction and partial death over several weeks, slowed mitochondrial movement in processes and appears to reproduce the Lewy neuritic changes of early Parkinson's disease pathology but does not cause Lewy body inclusions. The overall pattern of transcriptional regulation is gene under-expression with minimal regulation of ETC genes in spite of rotenone's being a complex I toxin. This rotenone-SH-SY5Y model in a differentiated human neural cell mimics changes of early Parkinson's disease and may be useful for screening therapeutics for neuroprotection in that disease stage.</p

Control of Canalization and Evolvability by Hsp90

Partial reduction of Hsp90 increases expression of morphological novelty in qualitative traits of Drosophila and Arabidopsis, but the extent to which the Hsp90 chaperone also controls smaller and more likely adaptive changes in natural quantitative traits has been unclear. To determine the effect of Hsp90 on quantitative trait variability we deconstructed genetic, stochastic and environmental components of variation in Drosophila wing and bristle traits of genetically matched flies, differing only by Hsp90 loss-of-function or wild-type alleles. Unexpectedly, Hsp90 buffering was remarkably specific to certain normally invariant and highly discrete quantitative traits. Like the qualitative trait phenotypes controlled by Hsp90, highly discrete quantitative traits such as scutellor and thoracic bristle number are threshold traits. When tested across genotypes sampled from a wild population or in laboratory strains, the sensitivity of these traits to many types of variation was coordinately controlled, while continuously variable bristle types and wing size, and critically invariant left-right wing asymmetry, remained relatively unaffected. Although increased environmental variation and developmental noise would impede many types of selection response, in replicate populations in which Hsp90 was specifically impaired, heritability and ‘extrinsic evolvability’, the expected response to selection, were also markedly increased. However, despite the overall buffering effect of Hsp90 on variation in populations, for any particular individual or genotype in which Hsp90 was impaired, the size and direction of its effects were unpredictable. The trait and genetic-background dependence of Hsp90 effects and its remarkable bias toward invariant or canalized traits support the idea that traits evolve independent and trait-specific mechanisms of canalization and evolvability through their evolution of non-linearity and thresholds. Highly non-linear responses would buffer variation in Hsp90-dependent signaling over a wide range, while over a narrow range of signaling near trait thresholds become more variable with increasing probability of triggering all-or-none developmental responses

The effects of acute tryptophan depletion on costly information sampling: impulsivity or aversive processing?

RATIONALE: The neurotransmitter serotonin (5-HT) has been implicated in both aversive processing and impulsivity. Reconciling these accounts, recent studies have demonstrated that 5-HT is important for punishment-induced behavioural inhibition. These studies focused on situations where actions lead directly to punishments. However, decision-making often involves making tradeoffs between small 'local' costs and larger 'global' losses. OBJECTIVE: We aimed to distinguish whether 5-HT promotes avoidance of local losses, global losses, or both, in contrast to an overall effect on reflection impulsivity. We further examined the influence of individual differences in sub-clinical depression, anxiety and impulsivity on global and local loss avoidance. METHODS: Healthy volunteers (N = 21) underwent an acute tryptophan depletion procedure in a double-blind, placebo-controlled crossover design. We measured global and local loss avoidance in a decision-making task where subjects could sample information at a small cost to avoid making incorrect decisions, which resulted in large losses. RESULTS: Tryptophan depletion removed the suppressive effects of small local costs on information sampling behaviour. Sub-clinical depressive symptoms produced effects on information sampling similar to (but independent from) those of tryptophan depletion. Dispositional anxiety was related to global loss avoidance. However, trait impulsivity was unrelated to information sampling. CONCLUSIONS: The current findings are consistent with recent theoretical work that characterises 5-HT as pruning a tree of potential decisions, eliminating options expected to lead to aversive outcomes. Our results extend this account by proposing that 5-HT promotes reflexive avoidance of relatively immediate aversive outcomes, potentially at the expense of more globally construed future losses

Re-analysis of the larval testis data on meiotic sex chromosome inactivation revealed evidence for tissue-specific gene expression related to the drosophila X chromosome

<p>Abstract</p> <p>Background</p> <p>Meiotic sex chromosome inactivation (MSCI) during spermatogenesis has been proposed as one of the evolutionary driving forces behind both the under-representation of male-biased genes on, and the gene movement out of, the X chromosome in <it>Drosophila</it>. However, the relevance of MSCI in shaping sex chromosome evolution is controversial. Here we examine two aspects of a recent study on testis gene expression (Mikhaylova and Nurminsky, <it>BMC Biol </it>2011, <b>9:</b>29) that failed to support the MSCI in <it>Drosophila</it>. First, Mikhaylova and Nurminsky found no differences between X-linked and autosomal genes based on the transcriptional profiling of the early testis development, and thus concluded that MSCI does not occur in <it>D. melanogaster</it>. Second, they also analyzed expression data from several <it>D. melanogaster </it>tissues and concluded that under-representation on the X chromosome is not an exclusive property of testis-biased genes, but instead, a general property of tissue-specific genes.</p> <p>Results</p> <p>By re-analyzing the Mikhaylova and Nurminsky's testis data and the expression data on several <it>D. melanogaster </it>tissues, we made two major findings that refuted their original claims. First, the developmental testis data has generally greater experimental error than conventional analyses, which reduced significantly the power to detect chromosomal differences in expression. Nevertheless, our re-analysis observed significantly lower expression of the X chromosome in the genomic transcriptomes of later development stages of the testis, which is consistent with the MSCI hypothesis. Second, tissue-specific genes are also in general enriched with genes more expressed in testes than in ovaries, that is testis-biased genes. By completely excluding from the analyses the testis-biased genes, which are known to be under-represented in the X, we found that all the other tissue-specific genes are randomly distributed between the X chromosome and the autosomes.</p> <p>Conclusions</p> <p>Our findings negate the original study of Mikhaylova and Nurminsky, which concluded a lack of MSCI and generalized the pattern of paucity in the X chromosome for tissue-specific genes in <it>Drosophila</it>. Therefore, MSCI and other selection-based models such as sexual antagonism, dosage compensation, and meiotic-drive continue to be viable models as driving forces shaping the genomic distribution of male-related genes in <it>Drosophila</it>.</p